Charged particle rapidity distributions at RHIC 
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Using a multiphase transport model (AMPT), whicli includes both initial partonic and final 
hadronic interactions, we study the rapidity distributions of charged particles such as protons, 
antiprotons, pions, and kaons in heavy ion collisions at RHIC. The theoretical results for the total 
charged particle multiplicity at midrapidity are consistent with those measured by the PHOBOS 
collaboration in central Au+Au collisions at = 56 and 130 AGeV. We find that these hadronic 
observables are much more sensitive to the hadronic interactions than the partonic interactions. 
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Collisions of nuclei at high energies offer the possibil- 
ity to subject nuclear matter to the extreme conditions 
of large compression and high excitation energies. Stud- 
ies based on both non-equilibrium transport models |^] 
and equilibrium thermal models have shown that the 
experimental data from heavy ion collisions at SIS, AGS 
and SPS, where the center of mass collision energies are, 
respectively, about 3, 5 and 17 AGeV, are consistent with 
the formation of a hot dense nuclear matter in the initial 
stage of collisions. With the Relativistic Heavy Ion Col- 
lider (RHIC) at Brookhaven National Laboratory, which 
can reach a center of mass energy of 200 AGeV, the ini- 
tial energy density is expected to exceed that for the 
transition from the hadronic matter to the quark-gluon 
plasma. Experiments at RHIC thus provide the oppor- 
tunity to recreate the matter which is believed to have 
existed during the first microsecond after the Big Bang 
and to study its properties. 

Recently, charged particle multiplicity near mid- 
rapidity has been measured in central Au-|-Au collisions 
at = 56 and 130 AGeV at RHIC by the PHOBOS col- 
laboration 1^ . The observed charged particle density per 
participant is found to be compatible with the predictions 
of the HIJING model that includes particle production 
from minijets produced in hard-scattering processes [Q. 
Although the HIJING model implements the parton en- 
ergy loss via jet quenching ||^, it does not include explicit 
interactions among minijet partons and the final-state in- 
teractions among hadrons. Other models have also been 
used to understand the data from the PHOBOS collabo- 
ration. The LEXUS model which is based on a linear 
extrapolation of ultra-relativistic nucleon-nucleon scat- 
tering to nucleus-nucleus collisions, predicts too many 
charged particles compared with the PHOBOS data [Q. 
On the other hand, the hadronic cascade model LU- 
CIFER j|] predicts a charged particle multiplicity near 
mid-rapidity that is comparable to the PHOBOS data 
1^. In this paper, we shall use a recently developed mul- 
tiphase transport model (AMPT) |]l^ , that includes both 
partonic and hadronic interactions, to study their effects 
not only on the total charged particle multiplicity but 
also on those of kaons, protons, and antiprotons. 



In the AMPT model, the initial conditions are ob- 
tained from the HIJING model Q by using a Woods- 
Saxon radial shape for the colliding nuclei and includ- 
ing the nuclear shadowing effect on parton production 
via the gluon recombination mechanism of Mueller-Qiu 
pl| . After the colliding nuclei pass through each other, 
the Gyulassy-Wang model ||l^ is then used to generate 
the initial space-time information of partons. In the de- 
fault HIJING, these minijet partons are allowed to lose 
energy via the gluon splitting mechanism and transfer 
their energies to the nearby strings associated with initial 
soft interactions. Such jet quenching is replaced in the 
AMPT model by explicitly taking into account parton- 
parton collisions via Zhang's Parton Cascade (ZPC) [p^ . 
At present, only gluon elastic scatterings are included, so 
the partons do not suffer any inelastic energy loss as they 
traverse the dense matter. After partons stop interact- 
ing, they combine with their parent strings and are then 
converted to hadrons using the Lund string fragmenta- 
tion model 0,11^ after an average proper formation time 
of 0.7 fm/c. Dynamics of the resulting hadronic matter is 
described by a relativistic transport model (ART) [ p^ , 
which has been improved to include baryon-antibaryon 
production from meson-meson interactions and their an- 
nihilation using cross sections given in Ref. [^,|8|. Also, 
the K* resonances are explicitly treated by including 
their production from pion-kaon and pion-rho scatterings 
p9| and the inverse reactions of decay and absorption. 

We first determine the parameters in the AMPT 
model by fitting the experimental data from central 
Pb-|-Pb collisions at center of mass energy of 17 AGeV 
po| . Specifically, to describe the measured net baryon 
rapidity distribution, we have included in the Lund 
string fragmentation model the popcorn mechanism 
for baryon-antibaryon production with equal probabili- 
ties for baryon-meson-antibaryon and baryon-antibaryon 
configurations jl^. Also, to account for pion and the 
enhanced kaon yields in the preliminary data from the 
same reaction, we have modified two other parameters in 
the Lund string fragmentation model, following the ex- 
pectation that the string tension is increased in the dense 
matter formed in the initial stage of heavy ion collisions. 
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In the Lund string fragmentation model as imple- 
mented in the JETSET/PYTHIA routine one first 
assumes that a string fragments into quark-antiquark 
pairs with a Gaussian distribution in transverse mo- 
mentum. Also, a suppression factor of 0.30 is used for 
strange quark-antiquark pair production relative to the 
light quark-antiquark pair production. Hadrons are then 
formed from these quarks and antiquarks. For a hadron 
with a given transverse momentum m± determined by 
those of its quarks, its longitudinal momentum is given 
by the splitting function , 



f{z) cx z~\l - z)''exp{-b m\lz). 



(1) 



where z is the light-cone momentum fraction of the pro- 
duced hadron with respect to that of the fragmenting 
string. Based on the Schwinger mechanism for particle 
production in a strong field, the production probability 
is proportional to exp(— 7rm^/«;), where k is the string 
tension, i.e., energy in a unit length of string. The aver- 
age squared transverse momentum of produced particle is 
thus proportional to k. Increasing the string tension then 
leads to a broader distribution of the transverse momenta 
of produced quark-antiquark pairs and also a reduced 
suppression for strange quark-antiquark pairs. Since the 
average squared transverse momenta of produced parti- 
cles obtained from Eq. ([^) is (p^) = [6(2-|-a)]^^ for mass- 
less particles, the two parameters a and h are approxi- 
mately related to the string tension by k (x [6(2 -I- a)]~^. 
In the HIJING model, the default values for a and h are 
0.5 and 0.9 GeV~^, respectively. We change their values 
to 2.2 and 0.5 GeV~^ in order to increase, respectively, 
the pion and kaon multiplicities. These values of a and 6 
correspond to a 7% increase of the string tension. 

Resuhs from the AMPT model for Pb+Pb coUisions 
at %fs = 17 AGeV, based on the above modification of 
the HIJING and ART models, are shown in Fig. ||. It is 
seen that our model gives a reasonable description of the 
data on the rapidity distributions of total and negatively 
charged particles, net-protons and antiprotons, charged 
pions, and charged kaons. Since the probability for mini- 
jet production is very small in collisions at SPS energies, 
the partonic stage does not play any role in the collisions. 
We find that final-state hadronic scatterings reduce the 
proton and antiproton yields, but increase the produc- 
tion of kaons and antikaons by about 20%. In contrast, 
kaon yields in the default HIJING model are smaller than 
our final results by about 40%. 

Since the number of strings associated with soft in- 
teractions in the HIJING model depends weakly on the 
colliding energy, the parameters in the string fragmenta- 
tion model are not expected to change much with the en- 
ergy. We thus use the same parameters determined from 
the experimental data at SPS to study heavy ion colli- 
sions at RHIC energies. In Fig. ^ the results for central 
Au-|-Au collisions at center of mass energies of 56 AGeV 
(dashed curves) and 130 AGeV (solid curves) are shown 
together with the data from the PHOBOS collaboration 



. The measured total charged particle multiplicities at 
mid-pseudorapidity at both energies are well reproduced 
by our model. 
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FIG. 1. Rapidity distributions of total and negatively 
charged particles (upper left panel) , net-protons and antipro- 
tons (upper right panel), charged pions (lower left panel), 
and charged kaons (lower right panel) in heavy ion collisions 
at = 17 AGeV. The circles are the experimental data for 
5% most central Pb-|-Pb collision from the NA49 Collabora- 
tion, and the solid curves are the AMPT model calculations 
for impact parameters of 6 < 3 fm. 
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FIG. 2. Same as Fig. ^ but for Au+Au collisions at 
^ = 56 and 130 AGeV. The solid circles are the PHOBOS 
data for 6% most central collisions while the curves are the 
AMPT calculations for fe < 3 fm. 
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The energy dependence of charged particle yields at 
mid-rapidity from the SPS to RHIC energies is shown 
in Fig. ^. The proton yield is seen to have a mini- 
mum at energies between SPS and the highest energy at 
RHIC, while antiproton yield increases almost linearly 
with Ins. As a result, the p/p ratio increases rapidly 
from about 0.1 at SPS to about 0.8 at the RHIC en- 
ergy of ^/s = 200 AGeV, indicating the formation of a 
nearly baryon-antibaryon symmetric matter at high en- 
ergies. Meson yields in general exhibit a faster increase 
with energy; in particular, we find that the K'^ ji:^ ra- 
tio is almost constant within this energy range, suggest- 
ing the approximate chemical equilibrium for strangeness 
production. The / ratio increases gradually from 
0.7 at SPS to about 1.0 at = 200 AGeV as a result of 
the nearly baryon-antibaryon symmetric matter formed 
at high energies. 
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FIG. 3. Energy dependence of charged particle yields at 
mid-rapidity. The ratios of K~ jK^ , p/p and /n^ are 
shown in the insert. 



To see the effects of hadronic interactions, we show 
in Fig. ^ by dashed curves the rapidity distributions of 
charged particles obtained from the AMPT model with- 
out the ART model for central Au-|-Au collisions at 130 
AGeV. In this case, there is a significant increase in the 
numbers of total charged particles, pions, protons, and 
antiprotons at midrapidity. The kaon number is, on the 



other hand, reduced slightly. As a result, the ratios of p/p 
and i<r+/7r+ in the absence of final-state hadronic interac- 
tions are 0.80 and 0.13, respectively, instead of 0.66 and 
0.18 from the default AMPT model. We note that al- 
though the defauh HIJING |] gives a total charged par- 
ticle multiplicity at midrapidity that is consistent with 
the PHOBOS data, including hadronic scatterings would 
reduce its prediction appreciably. 
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FIG. 4. Rapidity distributions of charged particles for 
central (6 < 3 fm) Au-I-Au collisions at 130 AGeV from 
AMPT model with default parameters (solid curves), without 
hadronic scatterings (dashed curves), and with jet quenching 
of dE/dx = 1 GeV/fm (dot-dashed curves). 

Effects of partonic dynamics on the final hadronic ob- 
servables can also be studied in the AMPT model. Turn- 
ing off the partonic cascade in the AMPT model, we 
find that this leads to less than ^ 5% change in the 
final charged particle yields at ^/s = 130 AGeV. This 
indicates that the multiplicity distribution of hadrons 
are not very sensitive to parton elastic scatterings. To 
take into account the effect of parton inelastic collisions, 
which are mostly responsible for energy loss, we include 
in the AMPT model also the default jet quenching, i.e., 
an energy loss of dE/dx = 1 GeV/fm, before minijet par- 
tons enter the ZPC parton cascade. The results with jet 
quenching for central Au-f Au collisions at 130 AGeV are 
shown in Fig. ^ by the dot-dashed curves. We see that 
the quenching effects are larger for pions than for kaons, 
protons, and antiprotons. Since the present calculations 
from the AMPT model without jet quenching already re- 
produce the data at the energy of 130 AGeV, and further 
inclusion of jet quenching of dE/dx = 1 GeV/fm entails 
a 10% increase of the final yield of total charged particles 
at midrapidity, our results for the rapidity distribution of 
charged particles are thus consistent with none or a weak 
jet quenching at this energy. 
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We note that without initial nuclear shadowing on par- 
ton production the charged particle multiplicity at mid- 
rapidity at 130 AGeV increases by about 30%. This in- 
crease can nevertheless be offset by using different val- 
ues for the parameters in the Lund string fragmenta- 
tion. Since nuclear shadowing has negligible effects at 
SPS energies due to insignificant production of minijets, 
to reproduce both SPS and RHIC data using the same 
parameters requires the inclusion of nuclear shadowing 
on parton production. 

In summary, using a multiphase transport model 
(AMPT), which includes both initial partonic and final 
hadronic interactions, we have studied the rapidity distri- 
butions of charged particles such as protons, antiprotons, 
pions, and kaons in heavy ion collisions at RHIC. With 
the model parameters constrained by central Pb-|-Pb col- 
lisions at ^/s = 17 AGeV at SPS, the theoretical results 
on the total charged particle multiplicity at midrapidity 
in central Au-|-Au collisions at yfs = 56 and 130 AGeV 
agree quite well with the data from the PHOBOS col- 
laboration. We find that the antiproton to proton ra- 
tio at mid-rapidity increases appreciably with ^/s, indi- 
cating the approach to a near baryon-antibaryon sym- 
metric matter in high energy collisions. Furthermore, 
the K'^ It:~^ ratio is almost constant within the energy 
range studied here, suggesting the approximate chemical 
equilibrium for strangeness production in these collisions. 
These hadronic observables are, however, less sensitive to 
the initial partonic interactions than the final hadronic 
interactions. To observe the effects of the partonic matter 
formed in the initial stage thus requires measurements of 
other observables such as J / ip suppression |2l[] , the ellip- 
tic flow ||2^ , and high p± spectra Q . The magnitude of 
elliptic flow has been shown to be sensitive to the parton- 
parton cross sections in the ZPC parton cascade model 
p4| , and the J/t/j suppression results using the AMPT 
model indicate that the partonic matter plays a much 
stronger role than the hadronic matter ||2^ . 
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